Magnetodielectric coupling of infrared phonons in single crystal Cu 2 OSe03 
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Reflection and transmission as a function of temperature have been measured on a single crystal 
of the magnetoelectric ferrimagnetic compound Cu2 0Se03 utilizing light spanning the far infrared 
to the visible portions of the electromagnetic spectrum. The complex dielectric function and optical 
properties were obtained via Kramers-Kronig analysis and by fits to a Drude-Lortentz model. The 
fits of the infrared phonons show a magnetodielectric effect near the transition temperature (T c ~ 
60 K). Assignments to strong far infrared phonon modes have been made, especially those exhibiting 
anomalous behavior around the transition temperature. 



I. INTRODUCTION 

The magnetodielectric effect refers to a change in the 
dielectric constant induced by an external magnetic field 
or by the onset of spontaneous magnetization. 1 . Observ- 
able anomalies in the polarization of a material near a 
magnetic transition indicate the presence of finite mag- 
netoelectric coupling, i.e., a cross coupling between mag- 
netic and electric ordersi 2 - Magnetoelectric coupling has 
attracted significant interest* 3 - - — recently not only for sci- 
entific purposes but also for its use in novel technolog- 
ical devices. One can envision field sensors and mag- 
netic memory switches where the magnetic order could 
be altered by adjusting the electric field. A shortcom- 
ing of magnetoelectrics has been the scarcity of materi- 
als that can simultaneously support magnetic and elec- 
tric orders. The fundamental issue behind this scarcity 
does not lie in the symmetry demands of the two orders; 
rather, the particular electron configurations that favor 
their origins are of a contradictory nature.— The elec- 
tron configurations required by magnetism usually pre- 
vail and the electric ordering is forced to originate from 
"improper" means. 7 Recent studies on materials such as 
TbMn 2 5 £ DyMn 2 5 £ and BiMnO^ all show magne- 
todielectric effects near magnetic transition temperatures 
where a distortion of the unit cell is also observed, as re- 
quired by magnetic ordering. This evidence suggests that 
the change in lattice constant, which results in a corre- 
sponding change in bond lengths, is the mechanism that 
alters the polarity at T c and, hence, produces a dielectric 
anomaly. 

Here we report our infrared studies on a single crystal 
of Cu20Se03, a piezoelectric with a ferrimagnetic tran- 
sition temperature of T c ~ 60 K. 1 Although we found no 
drastic anomalies across T c , a thorough inspection of the 
data combined with some modeling lead us to a mag- 
netodielectric effect. Recently Gnezdilov et al^ have 
presented a Raman study of Cu20Se03 prepared in the 
same way as our crystal. They observed the abrupt ap- 
pearance of 3 new lines in the spectra upon cooling below 
T C) and an additional 2 lines that appeared below 20 K. 
Gnezdilov et al. also gave a detailed description of the 
crystal and magnetic symmetry of this compound. 



II. EXPERIMENTAL PROCEDURES 

Single crystals of Cu2 0Se03 were grown by a stan- 
dard chemical vapor phase method. Mixtures of high pu- 
rity CuO (Alfa-Aesar, 99.995%) and Se0 2 (Alfa-Aesar, 
99.999%) powder in molar ratio 2:1 were sealed in the 
quartz tubes with electronic grade HC1 as the trans- 
port gas for the crystal growth. The ampoules were 
then placed horizontally into a tubular two-zone furnace 
and heated slowly by 50° C/h to 600° C. The optimum 
temperatures at the source and deposition zones for the 
growth of single crystals have been found to be 610° C 
and 500° C, respectively. After six weeks, many dark 
green, indeed almost black, Cu20Se03 crystals with a 
maximum size of 8 x 6 x 3 mm 3 were obtained. X- 
ray powder diffraction (XRD) analysis was conducted on 
a Rigaku X-Ray diffractometer with Cu Ka radiation 
(A = 1.5418 A). An electron microprobe was used for 
chemical analysis of all solid samples. 

Effenberger and Pertliki 2 - solved the crystal struc- 
ture using single-crystal X-ray diffraction. The compact 
crystal structure consists of three basic building blocks, 
square pyramidal CuOs, trigonal bipyramidal CuOs, and 
a lone pair containing tetrahedral Se03 unit. The oxy- 
gen atoms in the unit cell are shared amongst the three 
building blocks. The square pyramidal CuOs units ex- 
ist in a 3-to-l ratio to the trigonal bipyramidal CuOs 
units within the conventional unit cell. This ratio will 
be important subsequently when explaining the magnetic 
structure. All copper ions possess a +2 oxidation state. 
More detailed descriptions of the crystal structure are 
found elsewhereiii 3 - 

The material crystallizes in the P2i3 cubic space group 
and has been shown to remain metrically cubic with 
no abnormal change in the lattice constant through the 
Curie temperature and down to 10 K.- The onset of mag- 
netic order does have the effect of reducing the crystal 
symmetry to R3. Full cubic symmetry would require all 
copper ions to feel the same Coulomb interaction from 
nearest neighbor copper spins. The proceeding expla- 
nation is the case for ferromagnetism and antifcrromag- 
netism but not for ferrimagnetism, which is why a reduc- 
tion from cubic symmetry must accompany this ordering. 
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The temperature dependent (5-300 K) reflectance 
and transmittance measurements were collected using a 
Bruker 113v Fourier Transform interferometer in con- 
junction with a helium cooled silicon bolometer detec- 
tor in the spectral range 30-700 cm -1 and from 700- 
5,000 cm -1 with a nitrogen cooled MCT detector. Room 
temperature measurements from 5,000-40,000 cm -1 were 
obtained with a Zeiss microscope photometer. After mea- 
suring the bulk reflectance over the entire spectral range, 
the crystal was polished to a thickness of 194 fim for 
transmittance measurements. All measurements were 
performed using unpolarized light at near-normal inci- 
dence. Magnetic measurements were performed in a com- 
mercial SQUID magnetometer (Quantum Design MPMS- 
XL7) on a single crystal sample mounted with the [111] 
axis parallel to the applied field. After cooling the sam- 
ple in zero field to 50 K, magnetization was measured 
in an applied field of 10 G while warming to 70 K. The 
dc susceptibility was calculated in the low-field limit as 
x(T) = M(T)/H. In addition, the isothermal magneti- 
zation as a function of applied field was measured at a 
temperature of 2 K, while sweeping the field from zero 
to 2 kG and back to zero. 



III. RESULTS AND ANALYSIS 
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FIG. 1: The dc susceptibility, x(T), of Cu 2 OSe0 3 near the 
ordering temperature (T c = 60 K) in an applied field of 10 G 
as measured while warming after zero- field cooling to 50 K. 
The inset shows the isothermal magnetization, M(H), as a 
function of applied field at 2 K, where the field was swept 
up to 2 kG before being reduced to zero. In all instances, 
the lines connecting the data points are guides for the eyes. 
The schematic shows the orientation of the single crystal with 
respect to the field applied parallel to the [111] direction. 



A. Magnetism 



Recent studies have measured the magnetic suscepti- 
bility of powdered samples Cu20Se03, finding ordering 
temperatures of T™ &cction = 55 K 14 and T° nsot = 60 K ± 
Because anomalies in the infrared spectrum at the tran- 
sition temperature are important, an accurate determi- 
nation of T c for the single crystal of interest was desired. 
The measured dc susceptibility as a function of temper- 
ature, x{T), is shown in Fig. [TJ Taking the transition 
temperature to be where the susceptibility turns upward, 
yonset = 50 K is found, consistent with the observations 
of Bos et aL- At 2 K (Fig.Q]inset), well within the ordered 
state, the magnetization saturates in a field of 800 G at 
1.0 Nfis, which is half of the expected saturation value 
for aS= 1/2 system, indicating a ferrimagnetic order- 
ing in a three-up and one-down configuration. No coer- 
cive field was measured; however, an inflection point with 
some slight hysteresis was observed near 400 G, which is 
also consistent with the findings of Bos et a/.— 

B. Reflectance and Transmittance Spectra 

The temperature dependent reflectance spectrum of 
Cu 2 OSe0 3 between 30 and 1,000 cm' 1 (4-120 meV) is 
shown in Fig. [2J A strong sharpening of many phonon 
modes is observed with decreasing temperature. It 
should be noted that there are no drastic anomalies in 
the far-infrared spectrum, such as the presence of new 
modes or the splitting of existing modes, consistent with 



no change in lattice constant at T c . Infrared spectroscopy 
is extremely sensitive to changes in dipole moments, espe- 
cially those brought about by lattice distortions. There- 
fore, the absence of these anomalies gives strong sup- 
port to the assertion of no lattice distortions at T c as 
initially determined by X-ray diffraction measurements. 1 
The inset of Fig. [5] shows the 300 K reflectance up to 
40,000 cm . The low absorption in the infrared is con- 
sistent with the insulating nature of the compound. The 
onset of electronic absorption is indicated by the upturn 
of the spectra around 26,000 cm -1 (3.2 eV). 

The transmittance spectra, as depicted in Fig. El ar e in 
good agreement with the reflectance measurements. At 
frequencies below the strong phonon modes (< 80 cm -1 ), 
the crystal transmits. The periodic oscillations in this 
range are the Fabry-Perot interference fringes due to mul- 
tiple internal reflections. Transmission gaps open be- 
tween the infrared phonon modes. These regions be- 
come increasingly more evident as temperature is lowered 
and the modes sharpen. The low-frequency transmission 
spectrum exhibits a weak phonon with a resonance fre- 
quency of 68 cm -1 that first appears as a small structure 
around 120 K and strengthens with decreasing temper- 
ature. The mode has a small oscillator strength, thus 
explaining why it was not observed in reflection. The 
weak phonon shows no response to a field of 0.14 kG 
parallel to the crystal surface; however, a magnetic ori- 
gin for this structure cannot be ruled out and is worthy 
of future investigation. 
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FIG. 2: (Color online) Temperature-dependent reflectance 
spectrum of Cu2 0Se03. The inset shows the 300 K re- 
flectance out to 40,000 cm" . 



1,000 cm" 1 . The principal effect is a sharpening of most 
modes, yielding a larger conductivity at the resonant fre- 
quency. 
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FIG. 4: (Color online) Far Infrared optical conductivity 
Cu 2 OSe0 3 at 20 K(red) 70 K (black). 



The phonon modes in reflectance arc viewed as 
Lorentzian oscillators (derivative-like features) in the op- 
tical conductivity that makes them intuitive for model- 
ing as harmonic oscillators. The vanishingly small static 
limit of <Ti(w) and the low background level of conduc- 
tivity throughout the infrared regions is further evidence 
of the insulating nature of the compound. 
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FIG. 3: (Color online) Temperature dependent transmission 
spectrum of 194 fj,m thick Cu20Se03 single crystal below 
100 cm" 1 . This region is highlighted by a weak phonon (~ 
68 cm - ) that begins to appear around 120 K. The inset shows 
transmission at 300 K measured out to the 40,000 cm" . 



C. Kramers-Kronig analysis and Optical Properties 

Kramers-Kronig analysis can be used to estimate 
the real and imaginary parts of the dielectric constant 
from the bulk reflectance R(oj)r^ Before calculating the 
Kramers-Kronig integral, the low frequency data were 
extrapolated as a constant for u> — > as befits an insula- 
tor. At high frequencies the reflectance was assumed to 
be constant up to 1 x 10 7 cm" 1 after which R ~ (w)~ 4 
was assumed as the appropriate behavior for free car- 
rier. The optical properties were derived from the mea- 
sured reflectance and the Kramers-Kronig-derived phase 
shift on reflection; in particular, we estimated the real 
part of the optical conductivity, cti(cj). Figure 2] depicts 
cti(w) at 20 K (below T c ) and 70 K (above T c ) from 30 to 



D. Oscillator- model fits 

The Drude-Lorentz model was used to fit the re- 
flectance and obtain a second estimate of the complex 
dielectric function in the infrared range. The model con- 
sists of a damped oscillator for each putative phonon in 
the spectrum plus a high frequency permittivity Eqo that 
describes the contribution of all electronic excitations. 
The model has the following mathematical form: 



E 



ujj 2 — uj 2 — iujjj 



(1) 



where Sj, LJj, and jj represent the oscillator strength, 
center frequency, and linewidth of the jth damped oscil- 
lator. The complex dielectric function provided by the 
Drude-Lorentz model is used to calculate the reflectivity, 
which agrees well with our original measured quantity. 
Figure [3] shows the calculated reflectivity and the mea- 
sured reflectivity at 70 K. 



4 




200 400 600 

Frequency (cm -1 ) 



800 1000 



lowing simple formula for the static dielectric constant: 

oo 

£ ° = 3 i + £ °° ( 2 ) 

The calculated static dielectric constant at temperatures 
surrounding T Cl is shown in Fig. [5] There is an anomalous 
jump at around 60 K. It should be noted that this magne- 
todiclectric effect in Cu20Se03 was previously observed 
by Bos et alX through dielectric capacitance measure- 
ments. While the value of the dielectric constant mea- 
sured here is considerably smaller than that in the previ- 
ous report, the direction and magnitude of the anomaly 
at T c are in good agreement. 



FIG. 5: (Color online) Experimental reflectance (blue points) 
and calculated reflectance (red line) from the Drude-Lorentz 
model of the Cu20SeC>3 70 K dielectric function. 



O) 8.2 
Q 




40 60 80 

Temperature (K) 



FIG. 6: The static dielectric constant as calculated from the 
Drude-Lorentz model at temperatures between 20 and 100 K, 
including T c . 



B. Anomalous phonons 

The observed phonons in Cu20Se03 may be divided 
into two categories: conventional and anomalous. Con- 
ventional phonons show a slight hardening of their fre- 
quencies when cooled to low temperatures, accompanied 
by a significant reduction in linewidth and at most mod- 
est changes in oscillator strength. Moreover, the tem- 
perature variation is smooth with no sudden changes in 
slope or value. Anomalous phonons violate one or more 
of these expectations. We have identified 11 anomalous 
phonons in our spectra, which show 28 total modes. A 
numbering scheme has been used which corresponds to 
the sequence of appearance of oscillators in the infrared 
spectrum starting with #1 (lowest frequency mode) and 
ending with #28 (highest frequency mode) 

Figure [7] displays the temperature dependence of the 
three oscillator parameters for 3 of the 11 anomalous 
phonons in the infrared spectrum. The oscillator param- 
eters for one conventional phonon are shown for compar- 
ison. It was observed that anomalous behavior usually 
could be found in all three oscillator parameters. 



IV. DISCUSSION 

A. Magnetodielectric effect 

Equipped with oscillator parameters to describe each 
of the infrared phonons at all measured temperatures, 
one is now in position to closely monitor the subtle dy- 
namics of the phonon structures across the transition 
temperature. Despite the lack of drastic changes in 
the phonon spectrum at T c , it is worthwhile examining 
whether a combination of many small anomalies in the 
phonon dynamics might sum to give an overall effect. At 
this point it is logical to examine the static dielectric con- 
stant because it is the sum of parameters that describe 
the dielectric nature of the compound (the side to which 
the infrared is most sensitive) . Taking the zero frequency 
limit of the Drude-Lorentz formula, one arrives at the fol- 



C. Assignment of phonon modes 

Since the magnetodielectric effect is observed through 
lattice dynamics, analysis of the phonon spectrum is in- 
cumbent. The number of phonon modes expected in 
Cu20Se03 can be estimated by space group analysis. Us- 
ing the SMODES program^ we arrive at the following 
distribution of modes: 

Yopticai _ i^(R) _|_ i4£;*(- R ) -|- 4iy(^^-R) (3) 

where (R) and (IR) denote Raman active and infrared 
active modes respectively. We therefore have the poten- 
tial of 41 total infrared active modes, all of which pos- 
sess threefold degeneracy as indicated by their irreducible 
representation. However, only 28 modes in the infrared 
spectrum are detected, and a discussion about this dis- 
crepancy between the observed and predicted modes is 
beyond the scope of this work. 
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FIG. 7: (Color online) The (a) center frequency, (b) normalized oscillator strength, and (c) linewidth of four oscillators as a 
function of temperature. Oscillator #3 is a typical conventional phonon, whereas the other three oscillators show anomalous 
behavior as temperature is lowered across T c . 



We performed lattice dynamical calculations based on 
a real-space summation of screened coulomb interac- 
tions involving a spherical cut-off boundary^ Frequency, 
mode intensity, as well as displacement pattern were cal- 
culated based on the structure and valence as reported 
by Bos et alX Seven strong modes in the infrared spec- 
trum have been assigned by comparing the calculated 
and experimental spectra. 

Tabic U details the mode assignments made using the 
adopted numbering scheme and center frequency of each 
oscillator for identification. It should be noted that oscil- 
lators #12 and #15, which exhibit anomalous behavior 
at the transition temperature (#12 is shown in Fig. [7]), 
are associated with vibrations of oxygen around the cen- 
tral copper, the ion responsible for magnetic ordering. 



V. CONCLUSIONS 

Far-infrared measurements of single crystal reflectance 
from CU2OSCO3 reveal no drastic anomalies in the 
phonon spectrum as temperature traverses T c . However, 
a closer inspection of the dynamics of the phonon spec- 
trum, as modeled through Drude-Lorentz fitting, uncover 
an anomalous jump in the dielectric constant near T c . It 
is observed that 11 of the 28 total far- infrared phonons 



contribute to this magnetodielectric effect. A few strong 
far-infrared phonons have been assigned to motion of the 
CUO5 and Se03 units via a lattice dynamical calcula- 
tion. It is noteworthy that 2 of the 11 modes exhibiting 
anomalous behavior across T c have been assigned to the 
motions of oxygen around the central copper, the ion re- 
sponsible for magnetic ordering. A weak phonon that 
was not resolved in reflectance is observed below 120 K 
in the transmission spectrum. A magnetic origin for this 
structure has yet to be ruled out. Our infrared results 
agree with the Raman studies of Gnezdilov et alr^ in 
that we also observed a number of phonon modes that 
exhibit anomalies in their strength, center frequency, and 
linewidth, but we differ on other issues. For example, 
we do not observe any new modes below the magnetic 
transition whereas Genzdilov et ah 11 do detect some ad- 
ditional (rather weak and broad) features. They report 
three new modes appearing in the Raman spectra below 
the transition temperature at frequencies of ^261, 270, 
and 420 cm . The only structure we observe in these 
three spectral regions is at 270 cm -1 , where 1 of the 
11 reported anomalous phonons is present. Genzdilov et 
aZ^i also report two new lines originating below 20 K in 
the Raman spectra at ~86 and 203 cm -1 . Our infrared 
studies reveal a strong rather typical mode at ~86 cm -1 
and a weak anomalous mode at 203 cm -1 . If any new in- 
frared features possessed the same relative intensities as 
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TABLE I: Oscillator parameters (at 20 K) and their corresponding assignments for a few strong modes in the far infrared. 







20 K 




Cu 2 OSc0 3 




Index 


Osc Strength 


Center Freq 


Linewidth 


Anomalous 


Assignment 
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lo (cm -1 ) 


7 (cm" ) 






8 


0.521 


212 


2.9 


No 


be(J3 vibrating against CuOs 


12 


0.107 


336 


3.6 


Yes 


CuOs in plane vibration 


13 


0.188 


385 


2.0 


No 


SeC>3 bending mode 


15 


0.095 


437 


5.9 


Yes 


CuOrj out of plane vibration 


21 


0.291 


717 


15.8 


No 


SeC>3 antisymmetric stretch 


23 


0.054 


781 


13.9 


No 


Se03 antisymmetric stretch 


25 


0.016 


831 


5.1 


No 


Se03 radial breathing mode 


reported for the Raman spectra, we would have observed 
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